Lorentz Covariance and the Dimensional Crossover of 2d-Antiferromagnets 
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We derive a lattice /3-function for the 2d-Antiferromagnetic Heisenberg model, which allows the 
lattice interaction couplings of the nonperturbative Quantum Monte Carlo vacuum to be related 
directly to the zero-temperature fixed points of the nonlinear sigma model in the presence of strong 
interplanar and spin anisotropics. In addition to the usual renormalization of the gapful disordered 
state in the vicinity of the quantum critical point, we show that this leads to a chiral doubling of 
the spectra of excited states. 
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The 2d-Antiferromagnetic Heisenberg model is a well 
studied system that has two RG fixed points that have 
been identified at zero temperature throurii various ap- 
plications of the 2d nonlinear sigma model [l|, 0] . Depcn- 
dening on the relative anisotropy of the exchange cou- 
pling, J, between the two spatial directions of the system, 
the groundstate of the 2d antiferromagnetic groundstate 
is found to be either Neel ordered and gapless, or a gapful 
quantum disordered state. Varying the anisotropy drives 
dynamical fluctuations causing a zero ternperature phase 
transition between the two groundstates [3| , whereby the 
system effectively crosses over from 2d to 3d With- 
out the inclusion of a 0-term, the 2d nonlinear sigma 
model is used to give an effectively classical description 
of the Neel ordered groundstate, and the quantum na- 
ture of the disordered state arises in some nontrivial way 
through the dynamical scale evolution of the system. It 
has been questioned whether the inclusion of an explicit 
source term for quantum fluctuations would be of rele- 
vance in the Neel phase [1, [1] , in order to understand the 
mechanism of symmetry breaking, but the inclusion of 
such a 0-term provides an irrelevant perturbation. How- 
ever, if no source term for quantum fluctuations is in- 
cluded, the renormalization scale of the system can only 
be defined through phenomenological input. This pic- 
ture has been successfully verified in detail by comparing 
nonlinear sigma model predictions for the scaling of the 
correlation length in the Neel, and quantum disordered 
regimes, with numerics obtained from Quantum Monte 
Carlo (QMC) studies 0,11]. 

Recently a treatment of the reverse picture has been 
given via a conformal analysis of quantum spin chains 
[9|. Importantly, this identifies the effect of dynami- 
cal fluctuations on the 0-term. In principle, this effect 
should simply be to rescale the couplings of the under- 
lying sine-Gordon model by a dynamical factor. How- 
ever, we have recently pointed out that the convergence 
of the perturbative deformation is not guaranteed and a 
nonperturbative renormalization prescription is required 
(lo| . In this work, we extend this nonperturbative renor- 
malization picture to give a description of the role of 
quantum effects in the dynamical fluctuations of the 2d- 
Antiferromagnetic Heisenberg model. Whether or not a 



6'-term is relevant in the Neel phase [ll[, the 0(3) spin 
operators of the nonperturbative QMC method [I4I have 
already been obtained through a form of topological di- 
mensional reduction. Consequently, there are IR cutoff 
effects associated with the finite lattice system dynamics, 
and it is difficult to disentangle these effects from genuine 
mechanisms of symmetry breaking. Recent treatments 
based around the 2d 0(3) model have included the effect 
of a marginally irrelevant topological term possessing a 
U(l) symmetry Our new nonperturbative approach 
extends this picture by considering the effect of a finite 
IR lattice cutoff within the couplings of these terms. 

We focus on the nonperturbative properties of the 
continuous-time QMC method. This scheme has the 
same basic transfer matrix structure that is defined for 
the Density Matrix Renormalization Group (DMRG) 
method [l3|, [l3|, but the numerical basis of the lattice 
partition function for the continuous-time QMC method 
is defined through the application of classical loop-cluster 
methods. There are a number of closely related QMC 
schemes that can also be used to generate this lattice par- 
tition function, and these are discussed in detail in the 
review in [Tsj . However, we are single out the continuous- 
time scheme because the loop-cluster Monte Carlo updat- 
ing process is slightly different. The closed loops in this 
scheme, which represent the trace of the lattice parti- 
tion function, are generated through a sequence of suc- 
cessive local Monte Carlo updating decisions, and this is 
in contrast with other QMC schemes, closer to the origi- 
nal classical loop-cluster method, where the Monte Carlo 
updating decision is made by comparing the probabilistic 
weights of closed loops. All of the QMC schemes which 
use the DMRG form of the transfer matrix in [l^, [l3| 
can be given the same interpretation as systems of 0(3) 
spin vectors, but the continuous-time method has the ad- 
dtional feature that the lattice partition function is de- 
fined to be analytically continuous in Euclidcan-time over 
some finite interval. The reason why this is important is 
because, unlike the dynamical scaling relations of clas- 
sical spin systems [16| . the dynamical critical exponent 
for the QMC Euclidean-time direction only relates to a 
topologically complete system in the large lattice volume 
and stochastic probability distribution limits of the QMC 
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method [l7[. Therefore, a complete understanding of nu- 
merical correction effects and a direct comparison with 
the couplings of the 2d nonlinear sigma model is difficult 
with the generic method because finite-size lattice sys- 
tems only asymptotically approach these limits. Using 
the special continuous analytic property of the continous- 
time methods we are now able to define this limiting 
process in a new way. The nonperturbative lattice rcnor- 
malization group equation for the 2d-Antiferromagnetic 
Heisenberg model that, we now propose, is useful as it 
can save computational effort in numerical studies by in- 
terpolation of the dynamical scaling. 



I. QUANTUM MONTE CARLO METHOD 

In order to understand the renormalization group flow 
of the lattice couplings in the QMC method that are de- 
fined using the transfer matrices in [isl . fl^ . and how 
these relate to the 2d nonlinear sigma model, it is impor- 
tant to understand why these methods do not, in general, 
yield Lorentz invariant systems. This particular form of 
the QMC method has a deceptively simple overlap with 
the 2d 0(3) model. Numerical loop evolution generates 
the lattice path integrals and partition functions in the 
generic approach. This loop evolution proceeds over both 
the spatial and Euclidean-time extents of the lattice on 
an equal footing through the evaluation of local Monte 
Carlo decisions. The three component 0(3) quantum 
spin operators arc represented by their S'^-componcnt as 
a discrete spin defined on the spatial lattice sites, and the 
Sx and Sy components of the spin operators are defined 
through the projection of the probability distribution of 
the loop evolution defined for the Euclidean-time extent 
(ist . Thus we arrive at a very similar picture to the 
continuum 2d 0(3) model with a 0-term, but the cru- 
cial difference is that the topological term has a finite IR 
cutoff. Although for sufficiently large lattice volumes the 
distributions of the projections of the operators onto the 
spatial and Euclidean-time lattice directions do converge, 
these distributions are not constrained to be identical in 
the lattice ensemble through the definitions of the loop 
evolution (loj . Although the transfer matrix is isotropic 
in space and Euclidean-time, the individual loops that 
are realised through probabilistic decisions are not. Spa- 
tial and Euclidean-time isotropy is only realised in the 
stochastic limit of the lattice ensemble. In practice, a 
model defined with isotropic interaction couplings can 
be realised as being locally anisotropic through the dy- 
namics. 

Our aim in this article is to attempt to quantify the 
effect of this local anisotropy. In general, if we focus just 
on the spatial properties of the lattice ensemble then the 
critical behaviour, measured via the probability distri- 
bution of the S'z-component of spin, can be defined as 
a function of the lattice interaction coupling, J. This 
follows from the usual finite-size scaling (ESS) picture, 
which consists of asymptotic expansion about renormal- 



ization group fixed points [20| . However, since the generic 
QMC spin operators are also defined through a form of 
topological dimensional reduction, this implies that the 
fixed point can also be approached smoothly in (3 (the 
inverse temperature and cutoff scale for Euclidean-time). 
If this equivalence property is realised, the dynamical 
critical exponent is necessarily unity and the system is 
Lorentz invariant. However, in general, this scaling pic- 
ture is disrupted by IR cutoff effects, and consequently 
the numerically realised system cannot be Wick-rotated. 
Therefore by defining an exact basis of dynamical fluctu- 
ations, the equivalence with the 0(3) spin vectors is bro- 
ken in the generic QMC method that is defined through 
the transfer matrix in [l^, [3| . 

What we will now therefore do is simply change the em- 
phasis of the analysis. We will start with a Wick-rotated 
definition of the loop operators of the QMC method, to 
make the 0(3) equivalence exact, and then consider the 
IR cutoff effects that arise in finite-size lattice system in 
this choice of basis. This will enable us to quantify the 
breakdown of Lorentz invariance through the emergence 
of dynamical scale effects, and to quantify the effect of 
dynamical fluctuations on the irrelevant quantum pertur- 
bations of the 2d-Antiferromagnetic Heisenberg model 
through nonperturbative renormalization. We do this 



for the continuous-time QMC method in [13, because 
the lattice partition function is analytically continuous 
in Euclidean-time upto some finite scale. 



II. WICK ROTATION 

The new approach we take is to identify an analytic 
analogue of the QMC partition function, which is Wick- 
rotated. We express the 0(3) spin operators using a for- 
malism developed to treat the probabilistic dynamics of 
lattice systems analytically [2l|. This gives an exact de- 
scription of the lattice ensemble vacuum in the vicinity 
of the large volume stochastic limit, i.e. in the UV. We 
then expand this description towards the IR. The asymp- 
totic freedom of the 2d 0(3) model has been known for 
some time It is also known that the IR fixed point 
associated with the Berry phase can be treated via per- 
turbative deformations [2j|. What we treat here is the 
most general form of deformation of the irrelevant Berry 
term, via couplings which arc modified by the dynamical 
fiuctuations, where the convergence of the perturbative 
expansion would not then be guaranteed. 

In general, a nonperturbative lattice ensemble is only 
known through the expectation of its numerical matrix 
elements A. These matrix elements can be defined in 
terms of the generalised spin vectors n that describe the 
state of the spins on the lattice. For the QMC method, 
these matrix elements are defined in terms of both Sz 
components, defined on the spatial lattice sites, and also 
Sx and Sy components, defined through the projection of 
the probability distribution onto the Euclidean-time ex- 
tent of the lattice. In the stochastic limit we expect the 
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system to be Lorentz invariant, therefore we can define 
the Wick-rotation of the lattice system through the an- 
alytic continuation J = 16. The spin vectors defined on 
the spatial lattice sites will then represent the orientation 
of the Sx and Sy components of spin, and the spin vec- 
tors defined on the Euclidean-time direction will repre- 
sent the projection of the probability distribution defined 
for the Sz component of the spins. The reason for doing 
this is that in the continuous-time QMC method [l^l the 
number of lattice sites in the Euclidean-time direction is 
variable. Therefore, if we Wick-rotate the definitions we 
are able to consider the ESS of the topological terms and 
the IR cutoff effects in quantifiable lattice units. 

The transfer matrix of the 2d AFM is a 8 x 8 matrix 
[2^ . and so define the lattice matrix elements through 
three projection indices: onto the Euclidean-time direc- 
tion, and onto the two separate spatial directions. The 
implicit isotropy between the Sx and Sy components of 
spin defines a locally conserved current, 6, associated 
with each spatial lattice site. We therefore define two 
separate matrix elements for the two different spatial di- 
rections on the lattice. 



T(g)0i(g)e2 

(s,s',s") aeG 



{n(Slsa®ls'a®'i-s"c7\0l) 



(1) 



Instead of a single compact Berry phase term, in this 
nonperturbatively-motivated operator formalism we have 
two separate source terms in 6, one for each of the two 
anisotropic spatial directions. There is also an implicit 
cross-term in 9i ® 62, which comes from V . Erom this 
latter term we are able to generate local states like the 
plaquette ordered groundstates considered in [2^ because 
of the four-fold symmetry coming from the relative signs 
of 61 and ^2- However, we can also consider what hap- 
pens if we are unable to continue either 9i or 62 through 
2tt because of the finite IR cutoff. There will still be an 
equal number of instantons and anti-instantons within 
the vacuum, conserving the total topological charge, but 
if the rotational symmetry about the Sx and Sy spin com- 
ponent plane is lost, then the four- fold parity symmetry 
is broken. This symmetry is broken down to a two-fold 
symmetry such that the instantons arc no longer sym- 
metric under reflection about the Sz component of spin 
and their chirality is lost. In the 2d 0(3) model the Nccl 
phase is defined by having a larger correlation length than 
that of the dynamical fiuctuations. The same is true of 
the groundstate realised by the generic QMC method, 
but at finite Trotter number, the same is not necessar- 
ily true of correlations in the 47r rotational symmetry of 
the instantons away from the limit in which the lattice 
ensemble is Lorentz invariant. 

To simplify the cross-term of 6*1 002, we can define the 
projection of the 9i component of V , 



(s,s',s") creG ^ ' 



(2) 



where G is the discrete Z{2) algebra of the Sz spins, 
a is an element of G, 81, 02 = i , and T is the Trotter 
number of the Euclidean-time extent of the lattice. The 
partition function for the 2d-Antiferromagnetic Heisen- 
berg model is then given, in terms of these matrix ele- 
ments, as a path-integral over these two conserved cur- 
rents. 



Z= I d[9i]d[02] exp I A{n,)tei+B{n,)i92 

- Vin,) ds] (3) 



where V comprises the off-diagonal contributions to 
the spin operator basis, which is simply a general ma- 
trix element on the T 9i ® 82 lattice with no special 
symmetries, 



T^e^e, In® Is.® Is'. (B Is" a\n) 

{s,s',s") creG ^ ' ' 



C{n) = {n\V^{ns')\n) = ^ss' s"a{n){l.s' \n) 

(s,s\s") aeG 

T(giei®e2 

X ^ A,y,,"„(n)(n©l,^©l^v©l^^/-|n^.) (5) 

{s,s',s") creG 

Substituting, this then yields the partition function, 



(4) 



Z = y d[0i]d[02] cxp[ j A{ns)iei + ^Cins)i9i 

+ B{ns)ie2 - V'{ns) ds^ (6) 
where, = y(n,) - CK). 

III. LATTICE SPACING 

Our motivation for making a Wick-rotation is to quan- 
tify the effect of dynamical fluctuations on the irrele- 
vant Berry term in the 2d-Antiferromagnetic Hcisenberg 
model. In the above new Wick-rotated operator defini- 
tions, 9i and 92 are noncompact abelian variables, and 
the Euclidean-time extent is used to represent the multi- 
plicity of these phases. Consequently, each of the discrete 
intervals of the Euclidean-time extent now corresponds 
to a different 0-vacuum. To relate this new formalism to 
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the usual compact definition of topological charge (given 
for the 2d 0(3) model) we must somehow select one of 
these 6'- vacua over the others to be our reference compact 
sector. Conveniently, this choice of 0-vacuum naturally 
arises from the nonperturbative dynamics in the form 
of the numerical expectation value that is realised by 
the projection of the probability distribution defined for 
Euclidean-time [lo| . Stating this simply, each spin site on 
the 01 (g) O2 lattice has an associated matrix element 
A, and the Euclidean-time projection of A has a larger 
value in one of the discrete intervals of the Euclidean- 
time extent than in the others. This is, however, a spa- 
tial site-specific result. In principle, each spatial site can 
have the maxima of this projection arise in a different 
Euclidean-time interval. Different local sites can, there- 
fore, correspond to different 0-vacua, which presents a 
source of nonintegrablc singularity in analytically con- 
tinuing the lattice interaction coupling, J [2^. These 
singularities, though, are precisely the form of dynamical 
fluctuation-induced IR cutoff effect that we are aiming to 
quantify. 

To quantify these properties, we introduce local mea- 
sures of lattice spacing, a and &, defined in units of 6 
and /3, respectively. These are given as the difference be- 
tween the projection of the matrix element realised on a 
given lattice site and the expectation of the projection 
averaged over the corresponding lattice extent, 

6ias' = A{ns') - {A{ns,))e^ 

(3 a, = A{n,) - {A{n,))T 
02 bs" = B{ns») - (B{n,„))e, 

Pb, = Bins) - {B{n,))T 

(7) 

Practically, A{ns) and B{ns) can be found directly 
as they correspond to the diagonal entries of the nu- 
merical transfer matrix of the continuous-time QMC 
method. Similarly, the two operator projections A{ns') 
and B{nsii) can be calculated by summing the local 
continuous- time QMC transfer matrix entries over the 
spatial sites rather than Euclidean-time. Thus, in a for- 
mal mathematical sense these lattice spacing definitions 
define the support of the matrix elements found by in- 
tegrating the local T 81 02 matrix elements over a 
given lattice direction. 

A related quantity is the parallel transport between 
neighbouring lattice sites for which we can define the 
following local derivative operators, 

dsA{ns) = A{ns+i) - A{ns) = iOiAins) 
dsB{ns) = B{ns+i) - B{ns) = i02B{ns) (8) 

IV. LATTICE /3-FUNCTION 

To identify a new nonperturbative lattice /^-function 
we define an effective action for the partition function 



of (6), which is given as an expansion in the self-energy 
terms of the dynamical basis. A general loop generated 
by the Monte Carlo process is of the form T (g) 0i (g) 02 
and is given by the action of y(n)-operators on the vac- 
uum. Similarly, the action of the yl(n)-operators gener- 
ates loops of the form T ® 0i, and the _B(n)-operators 
generate loops of the form T ©2. Thus, the self-energy 
terms of the Wick-rotated dynamical basis we have de- 
fined, which describe the projection of the A and B oper- 
ators onto V , are of the form 0i and 02, and correspond 
to the 0-symmetry breaking component of the vacuum. 
The effective action is of the form, 

S = Si + S2 + Si2 (9) 

where, 

51 = ds dsA{ns) - ds[V'ins)A{ns)] - V'{ns) 

Jo 

52 = / ds dsB{ns) - ds[V'{ns)B{ns)] - ^'(n,,) 

Jq 

^12= / ds dsVCius) - ds[V'ins)VC{ns)] + l/'(n,) 
Jo 

(10) 

The terms of the effective action are of the same form 
as the undcformcd action in ([3]). These terms are de- 
fined through a form of topological dimensional reduc- 
tion where an explicit integration is performed over P 
to project out the Euclidean-time dependence. Next we 
integrate out the effective action terms in a second topo- 
logical dimensional reduction step by integrating out the 
expansion over the compact sectors that define each lat- 
tice matrix element. The purpose of this second step is 
to make the effective action compact, which allows us to 
then make a direct comparison between our new RG flow 
and that of the 2d nonlinear sigma model. To do this we 
introduce two new variables; x which is a general posi- 
tion index on T 0i, and y which is a general position 
index on T 02 . 

Our aim, following Refi, is to treat the effect of ir- 
relevant quantum fluctuations on the relevant dynami- 
cal fluctuations of the 2d nonlinear sigma model, given 
m RefQ. Our new operator formalism is defined to be 
exact in the basis of these dynamical fiuctuations from 
the properties of the continuous-time QMC method. In 
practice, the numerics suffer from IR cutoff effects but 
we have Wick-rotated the operator definitions in order 
to quantify this effect on the topology of the 0(3) spin 
vectors. In principle, the IR cutoff implies that the 27r- 
rotational symmetry of the Sx and Sy plane components 
is broken. However, we have the freedom, from the renor- 
malization group properties of the relevant dynamical 
fluctuations, simply to rescale the couplings of the op- 
erators such that the 27r-rotational symmetry is present. 
There is then a new singularity that sits on the lat- 
tice vertices, which appears essentially by modifying a 
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finite number of tlie poles of the lattice system to branch 
points. 

The second topological dimensional reduction step is 
straightforward to evaluate for the first two terms in (9) 
using the local lattice spacing definitions given in (7). 

Si ^ (fx as' dlA{n) - as dlV {n) - as /\ Qs' 

Jo J~TT 

(11) 



52 = / / d'y hs" d'yBin) ~ bs d^V'{n) - 6, A 6,. 

Jo J-TT 



(12) 



gap above the ordered groundstate. What we should 
therefore find is that the spectrum of excited states ex- 
hibits a doubling due to the degeneracy of the broken 
chirality of the instanton anti-instanton pairs, when we 
are in the region of the classical dimensional crossover 
of the 2d quantum antiferromagnet with the presence of 
strong anisotropics due to these quantum fluctuation ef- 
fects. 

The sine-Gordon renormalization group equations for 
pT|) and (fT^l) are then trivially modified in the presence of 
anisotropy [28[. However, this leads to the emergence of 
a new unstable fixed point and finite renormalized region 
in the vicinity of the second order quantum critical point. 



These isotropic terms are of the form of the corresponding 
sine-Gordon model description of quantum spin chains. 
This is expected since the interaction between the two 
spatial lattice directions is ignored. These terms imply, 
following RefH, that the effect of quantum fluctuations 
can be simply rescaled the into the existing couplings, 
which describe the dynamical fluctuations. The cross- 
term, which links the two spatial lattice directions, is 
slightly more involved, however, since it involves a choice 
of branch. 

Sv2 = f ! (fx - [1 - (v/C(^)eJ dlV'in) 

+ {l + bsA bs'>)[l - (F'(n,))eJ dliy^] 
-9,(v/CK))ei5,(F'K))e, 

= f d?x as- h{l + bsAbs")dl[^C{n)] 

Jo J-TT 

-as A {bs A bs")dlV'{n) - (a, A a,-) A {bs A bs") 

(13) 



das 2 f - bs A bs" \ ^, 
^ = "2"^(, j 

^ = {as' -bshbs"){2-2as) 

The axes-crossing in the usual hyperboloid scaling pic- 
ture of the {a ST as') phase plane is simply shifted by the 
inclusion of branch singularities. This leads to a very 
simple analytic rescaling of the relevant couplings of the 
gap state, but it is one which is not easily quantified 
perturbatively in terms of a global cross-term (as in ll|) 
because of the local site dependence of the branch ]VJ\ . 
It was argued in [13] that such a contribution should be 
responsible for the appearance of deconfined spinous at a 
finite energy gap above the ordered groundstate. What 
we should, therefore, find is that the spectrum of ex- 
cited states exhibits a doubling due to the degeneracy of 
the broken chirality of the instanton/anti-instanton pairs 
when we are in the region of the classical dimensional 
crossover of the 2d quantum antiferromagnet through 
the Lorentz covariance of the irrelevant quantum fluc- 
tuations. 



The cross-term contribution to the effective action, 
within the contour between tt and — tt, is of the same form 
as (|1H) . but the contribution from the branch point on 
the boundary leads to an additional phase contribution 
to the action [13]. By comparing (fTTj) and p3)) we find 
that the cross-term contribution is implicitly anisotropic 
and, therefore is of the form of the doubel Sine-Gordon 
model not the sine-Gordon model 0]. This means that 
although the effect of quantum fluctuations on the dy- 
namical fluctuations is irrelevant (amounting simply to a 
change in the renormalization scale of the vacuum), the 
effect of quantum fluctuations is relevant for the Lorentz 
covariance scale of the dynamical fluctuations. The ef- 
fect of having a branch point associated with parity in 
the cross-term, is that the parity of the vacuum (that 
allows an exact symmetry between the instantons and 
anti-instantons) can be explicitly broken via quantum 
fluctuations. It is only in the limit when bs A bs" = 1 
that the effective action term in (jl3p is Lorentz invariant 
and this effect is vanishing. 



V. SUMMARY 

From early spinwave analyses of the 2d quantum an- 
tiferromagnet it was concluded that the role of quantum 
fluctuations is irrelevant to the stability of the Neel or- 
dered groundstate [1^. Subsequent refinement of these 
arguments has suggested that the crucial analytic prop- 
erty of the analysis, which enables this stability, is the 
continuity of generalised spin operators that describe the 
vacuum [5|, |30[. In this article, we have now considered 
the role of dynamical cutoff effects on these quantum 
ffuctuations in an exact basis of dynamical fluctuations 
following the nonlinear sigma model treatment given in 
Refill. The crucial difference from previous studies is 
that we have now considered the role of anisotropy on 
modifying the renormalization cutoff scale on the phe- 
nomcnological couplings of the ordered state. Whilst it 
is known to high accuracy that the Neel ordered ground- 
state closely follows spinwave theory predictions [sij . 



6 



considerably less is known about relevant renormaliza- 
tion of the dynamical couplings by quantum fluctuations 
when the Lorentz symmetry of the system is broken. This 
nonperturbative scaling treatment forms a more realistic 
picture of the numerical scaling of the dynamical basis 
defined by the generic QMC method, which is only truly 
Lorentz invariant asymptotically close to the stochastic 
limit. Similarly, the new nonperturbative renormaliza- 
tion group formalism provides a loop expansion formal- 
ism which is suitable for rigorously probing the conformal 
correspondence of experimental high temperature super- 
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